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COMPUTER  MODELS  FOR  RAINFALL-RUNOFF 
AND  RIVER  HYDRAULIC  ANALYSIS1 
by 

Darryl  W.  Davis2 

INTRODUCTION 

The  application  of  computer  technology  to  analysis  of  the  rainfall- 
runoff  process  and  the  hydraulics  of  natural  rivers  has  greatly  expanded 
in  the  past  few  years.  A  large  number  of  special  purpose  programs  and 
a  few  programs  designed  for  general  application  have  been  developed  and 
applied  to  hydrologic  engineering  problems.  The  Hydrologic  Engineering 
Center  (HEC)  has  developed,  over  the  past  8  years,  a  number  of  generalized 
computer  programs  for  use  by  the  US  Army  Corps  of  Engineers  in  analyzing 
hydrologic  engineering  problems.  This  paper  briefly  describes  the 
activities  of  The  Hydrologic  Engineering  Center  and  discusses  the  capa¬ 
bilities  of  two  of  these  programs:  (1)  Flood  Hydrograph  Package  (HEC-1) 
and  (2)  Water  Surface  Profiles  (HEC-2) . 

THE  HYDROLOGIC  ENGINEERING  CENTER 

The  Hydrologic  Engineering  Center,  established  in  1964,  has  three 
principal  missions:  (a)  conduct  research  and  development  of  hydrologic 
engineering  techniques  for  use  in  the  Corps'  day-to-day  work,  (b)  provide 
training  for  Corps  of  Engineers  employees  in  traditional  as  well  as  newly 

^or  presentation  at  the  Spring  1973  Conference  of  the  Civil  Engineering 
Program  Applications  (CEPA)  Conference,  San  Francisco,  California 
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developed  hydrologic  engineering  techniques,  and  (c)  provide  assistance 
to  Corps  field  offices  in  hydrologic  engineering  studies.  The  Center 
has  recently  been  assigned  responsibility  in  the  area  of  planning 
analysis  so  that  its  capabilities  in  the  field  of  systems  analysis  can  be 
focused  on  priority  components  of  planning  problems.  From  the  beginning, 
HEC  concentrated  on  computer  applications  and  utilization  in  carrying 
out  each  of  these  missions  (reference  1). 

The  research  and  development  work  has  resulted  in  the  coding,  testing 
and  documentation  of  about  30  generalized  computer  programs  in  hydrologic 
engineering.  The  programs  are  generalized  in  that  they  can  be  applied 
without  modification  to  almost  any  problem  in  their  specific  area  of  appli¬ 
cation  regardless  of  the  scope  of  the  problem,  the  geographic  location  of 
the  problem,  or  the  degree  of  detail  required  in  a  particular  solution. 

Five  of  the  programs  are  large  programs  that  combine  a  number  of  functions 
into  a  single  package.  These  five  programs  are  in  the  areas  of  flood 
hydrograph  analysis,  water  surface  profiles,  reservoir  systems  analysis 
for  flood  control  and  conservation  and  monthly  streamflow  simulation. 

The  training  activities  conducted  or  sponsored  by  the  Center  are 
designed  to  Improve  the  hydrologic  engineering  capabilities  of  the  Corps 
of  Engineers  in  accomplishing  its  civil  works  mission.  This  specialized 
training  contributes  to  more  efficient  performance  of  technical  studies 
associated  with  the  planning,  design  and  operation  of  civil  works  projects. 
The  training  program  includes  about  20  weeks  per  year  of  1-  and  2-veek 
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courses  with  classroom  instruction  in  topics  of  hydrologic  engineering. 
Individual  training  and  short  seminars  are  also  included  when  desirable. 

The  special  assistance  mission  of  the  Center  serves  to  point  up  the 
needs  for  training  as  well  as  for  research  and  provides  realistic  problems 
for  testing  techniques  and  programs. 

THE  FLOOD  HYDROGRAPH  PACKAGE  (HEC-1) 

General  Capabilities 

HEC-1  can  be  characterized  as  a  single  storm  event  flood  runoff 
simulation  model.  Most  ordinary  flood  hydrograph  computations  associated 
with  precipitation  and  runoff  on  a  complex  multisubbasin,  multichannel 
river  basin  can  be  accomplished  with  the  program.  Because  of  the  modelling 
capability  of  the  program,  a  number  of  other  routines  ere  included  that  can 
assist  in  determining  the  appropriate  parameters  needed  to  model  the  runoff 
process  and  to  evaluate  the  effect  of  management  alternatives. 

The  five  major  types  of  computations  that  can  be  performed  by  HEC-1 

are: 

s  Generalized  precipitation,  runoff,  routing  and  combining 
operations  to  simulate  the  hydrologic  response  of  a 
watershed  and  its  stream  network  (the  modelling  element); 

•  Optimization  of  routing  parameters  (assistance  in  parameter 
derivation) ; 

•  Optimization  of  unit  hydrograph  and  loss  rate  parameters 
(assistance  in  parameter  derivation); 
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•  Scream  system  computations  for  specified  precipitation 
depth-area  storm  relationships  for  the  entire  watershed  or 
region  (application  of  modelling  capability) ; 

•  Specialized  precipitation  streamflow  network  simulation 
relative  to  multiple  floods  for  multiple  plans  of  basin 
development  and  the  economic  analysis  of  flood  damages 
(evaluation  of  management  alternatives) . 

In  the  process  of  modelling  a  basin,  the  program  provides  several 
techniques  for  inputting  and  distributing  the  precipitation,  treating 
the  precipitation  as  rainfall  or  snowfall,  computing  rainfall  and  snow¬ 
melt  losses  and  excess,  determining  subbasin  outflow  hydrographs  from  unit 
graph  techniques,  and  routing  hydrographs  by  storage  routing  methods.  If 
necessary,  different  techniques  for  each  process  may  be  combined  in  the 
same  job  for  the  basin  being  modelled.  Graphical  display  of  intermediate 
or  summary  hydrographs  and  precipitation  can  be  produced  where  desired. 

The  program  may  be  used  to  optimize  specified  parameters  of  the  precipi¬ 
tation  runoff  or  routing  processes  for  a  stream  reach  or  subbasin  to  achieve 
a  best-fit  with  respect  to  an  observed  hydrograph  and  known  precipitation 
or  a  known  inflow  hydrograph. 

The  precipitation  depth-area  stream  system  computation  is  designed  to 
compute  a  consistent  hydrograph  at  all  desired  points  in  a  complex  basin  so 
that  each  corresponds  to  a  specific  depth-area  relationship.  The  procedure 
operates  by  computing  simultaneously  a  maximum  of  five  base  floods,  each 


representing  average  rainfall  intensity  corresponding  to  a  specified 
area  size.  An  interpolated  hydrograph  is  automatically  established  for 
each  concentration  point  based  on  the  size  of  the  area  tributary  to  that 
point.  This  routine  is  useful  in  stream  systems  or  in  urban  storm 
drainage  computations  where  it  is  desired  to  compute  a  number  of  synthetic 
events  feuch  as  a  100-year  flood)  at  a  large  number  of  locations. 

The  routine  for  evaluating  reservoir  and  channel  development  plans 
for  one  or  more  locations  includes  the  computation  of  average  annual 
dollar  benefits  at  each  damage  center  for  each  plan  of  development  as  well 
as  for  existing  conditions.  This  involves  simultaneously  computing  a  nuober 
of  system  floods  for  each  plan,  covering  the  entire  range  of  floods  that 
significantly  contribute  to  damages.  The  floods  may  be  either  multiples  of 
the  runoff  from  a  single  representative  storm  or  the  runoff  from  multiples 
of  a  typical  storm  rainfall  pattern.  Flow-damage  relations  for  each  type 
of  damage  and  flood-peak  frequency  relations  for  existing  conditions  must 
be  specified  for  each  damage  center.  Unit  hydrograph  coefficients,  loss 
coefficients,  degree  of  imperviousness  and  routing  coefficients  for  each 
plan  must  also  be  specified. 

Modelling  Rainfall-Runoff  with  HEC-1 
The  program  is  designed  to  simulate  the  storm  rainfall-runoff 
process  and  is  composed  of  the  appropriate  mathematical  relationships  and 
conatants  that  describe  the  response  of  the  watershed  to  rainfall.  The 
model  accepts  total  storm  rainfall  for  each  subbasin,  deducts  losses  to 
determine  rainfall  excess,  transforms  the  excess  to  streamflow  by  the 
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unit  hydrograph  technique,  and  translates  the  streamflow  through  valley 
storage  and  reservoir  storage  by  simple  routing  procedures. 

Modelling  a  complex  basin  consists  of  defining  Its  topologlc  structure 
(subbasin  boundaries  and  areas,  stream  channels  and  the  logical  relation¬ 
ships  between  the  subbasins  and  stream  channels)  ,  and  defining  the  parameters 
that  describe  the  rain fall- runoff  response  of  the  subbasins  and  channels 
of  the  river  basin.  Parameters  are  needed  to  determine:  (1)  basin  average 
total  precipitation  and  its  time  distribution  for  each  subbasin,  (2)  precipi¬ 
tation  loss  rates  for  selected  storm  events,  (3)  unit  hydrographs  for  each 
subbasin,  (4)  base  flow  controls  for  each  subbasin,  and  (5)  routing  criteria 
for  each  channel  reach.  The  detailed  computational  algorithms  can  be  found 
in  reference  2.  Figure  1  is  a  simplified  flow  chart  of  HEC-1. 

Topology .  The  topology  of  the  basin  is  described  in  the  program 
essentially  by  the  way  the  sequences  of  operations  are  specified.  Hydrographs 
are  computed,  routed  and  combined  with  other  hydrographs  in  accordance  with 
the  data  sequence  provided.  In  this  manner  any  complex  basin  consisting  of 
a  very  large  number  of  subbasins  and  routing  reaches  can  be  accurately  described. 

Precipitation.  The  precipitation  applied  to  a  subarea  for  runoff  compu¬ 
tations  may  be  determined  by  three  methods.  The  computed  precipitation  will 
then  be  treated  as  rainfall  or  snowfall  depending  upon  temperature  criteria. 

a.  Subbasin  total  precipitation  is  computed  from  nonrecording 
station  precipitation  according  to  weights  provided  for  each  station.  The 
subbasin  total  precipitation  can  also  be  computed  considering  longer  term 
weighting  values  such  as  normal  annual  amounts  that  may  reflect  consistent 


Figure  I 

HEC-I  MACRO-FLOW  CHART 


•  Numbers  In  parentheaea  Indlcsta  the  value  of  firo*ra*  variable  JgPtP.  which  la  ured 
Co  describe  the  flva  bealc  typaa  of  Jobs  aa  follow* : 

(IJ  Routing  Ope  Imitation 

(2)  Loaa  Rat*  and  Unit  Graph  Optimisation 

())  Canaraliiau  Stream flow  ..acwork  Computations 

(4)  uapth-Ara*  Storn  Stream  System  Computation* 

(SI  feltlpla  Plan,  iultlpla  Ratio  Economic  Computation* 

♦  Miabora  In  bracket*  (  }  Indicate  tha  value  of  variable  IClHF  on  tha  K-catd  which  la  tha  kay 
aodal  building  card  describing  tha  currant  operation  to  be  undertaken. 

[0]  compute  runoff  froa  a  subarea 
ill  route  hydrograph  through  a  river  reach 
(2-*;j  combine  ICRHP  hydrograpba  at  thla  atatlan 

l Ml  compare  computed  and  observed  Hydrograph*  at  thla  atatlan 
(99]  «nd  or  job— print  summary  Information 


8 


/ 


physiographic  effects  in  the  subbasin.  The  subbasin-mean  recording 
precipitation  pattern  (time  distribution)  is  computed  in  a  similar 
manner  with  station  weights  as  appropriate.  The  subbasin-mean  precipita¬ 
tion  distribution  (time  distribution  of  total  storm  rainfall)  is  then 
computed  using  the  relative  pattern  to  distribute  the  total  subbasin 
rainfall. 

b.  Known  temporal  and  spatial  precipitation  for  a  subbasin  can 
be  supplied  directly.  By  this  method  precipitation  can  be  supplied 

for  each  Interval  or  as  a  time  pattern  for  a  given  storm  total.  Precipita¬ 
tion  distributions  may  also  be  specified  for  two  differently  sized  areas  and 
the  program  will  perform  a  logarithmic  interpolation  of  the  two  patterns 
for  a  specified  area. 

c.  The  program  provides  for  automatic  computation  of  standard 
project  storm  (SPS)  precipitation  using  Corps  of  Engineers  criteria  and 
for  probable  maximum  precipitation  using  criteria  developed  by  the  National 
Weather  Service.  For  the  SPS,  the  largest  day  of  precipitation  is  preceded 
by  the  second  largest  and  followed  by  the  third  largest.  Six-hour  storm 
amounts  within  each  day  are  similarly  distributed.  A  storm  transposition 
coefficient  can  be  supplied  or  will  be  computed  by  the  program  as  a  default 
option. 

Where  snowfall  and  snowmelt  are  considered,  there  is  provision  for 
separate  computation  in  up  to  ten  elevation  zones.  These  zones  sre  usually 
considered  to  be  in  elevation  increments  of  1,000  feet,  but  any  equal 
increments  of  elevation  can  be  used  as  long  as  the  air  temperature  lapse 
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rate  corresponds  to  the  change  in  elevation  within  the  zone.  The  input 
temperature  data  are  those  corresponding  to  the  bottom  of  the  lowest 
elevation  zone.  Temperatures  are  reduced  by  the  lapse  rate  in  degrees 
per  increment  of  elevation  zone.  The  base  temperature  at  which  melt  will 
occur  must  be  specified  because  variations  from  32°F  (0°C)  might  be 
warranted  considering  both  spatial  and  temporal  fluctuations  of  temperature 
within  the  zone. 

Precipitation  is  assumed  to  fall  as  snow  if  the  zone  temperature  is 
less  than  the  base  temperature  plus  2  degrees.  Melt  occurs  when  the 
temperature  is  equal  to  or  greater  than  the  baBe  temperature.  Snowmelt 
is  subtracted  and  snowfall  is  added  to  the  snowpack  in  each  zone.  Snow¬ 
melt  may  be  computed  by  the  degree-day  or  energy-budget  methods. 

Precipitation  Loss  Rates.  Loss  rates  can  be  computed  using  initial 
and  uniform  losses  or  by  a  function  which  relates  loss  rates  to  rainfall 
and  snowmelt  intensity  and  to  accumulated  loss  (ground  wetness) .  Figure  2 
shows  the  loss  rate  function  for  a  snow-free  basin.  The  loss  rate  function 
is  successively  applied  for  each  computational  Interval. 

The  loss  rate  parameters  of  figure  2  are: 

DLTKR  -  Amount  of  initial  accumulated  rain  loss  during  which 
the  loss  rate  coefficient  is  increased.  This  param¬ 
eter  is  considered  to  be  a  function  primarily  of 
antecedent  soil  moiBture  deficiency  and  is  usually 
different  for  different  storms. 
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STRKR  -  Starting  value  of  loss  coefficient  on  exponential  reces¬ 
sion  curve  for  rain  losses  (snow-free  ground) .  The 


starting  value  is  considered  a  function  of  infiltration 
capacity  and  thus  depends  on  such  basin  characteristics 
as  soil  type,  land  use  and  vegetal  cover. 

RT10L  -  Ratio  of  rain  loss  coefficient  on  exponential  loss  curve 
to  that  corresponding  to  10  inches  more  of  accumulated 
loss.  This  variable  may  be  considered  a  function  of  the 
ability  of  the  surface  of  a  basin  to  absorb  precipitation 
and  should  be  reasonably  constant  for  large  rather 
homogeneous  areas. 

ERAIN  -  Exponent  of  precipitation  for  rain  loss  function 

ALOSS  -  (AK  +  DLTK)  PRCPERAIN 

that  reflects  the  influence  of  precipitation  rate  on 
basin-average  loss  characteristics.  It  reflects  the 
manner  in  which  storms  occur  within  an  area  and  may  be 
considered  a  characteristic  of  a  particular  region. 

Varies  from  0.0  to  1.0.  The  terms  in  the  equation  are 
defined  as: 

ALOSS  «  loss  rate  for  particular  time  interval 
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AK 
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DLTK 


■  loss  rate  coefficient  at  beginning  of  time 
Interval,  value  on  STRKR  exponential  loss 
curve. 

rainfall  intensity  in  inches  (mm)  per  hour, 
incremental  increase  in  loss  rate  coefficient. 
DLTK  is  assumed  to  be  a  parabolic  function  of 
the  accumulated  loss  for  DLTKR  amount  of 
accumulated  loss.  DLTK  is  a  maximum  of  0.2 


DLTKR  initially  reducing  to  zero  when  the 
accumulated  loss  equals  DLTKR. 


Unit  Hydrograph.  The  unit  hydrograph  corresponding  to  the 
appropriate  duration  of  rainfall  excess  can  be  supplied  directly  or  nay 
be  conputed  from  coefficients  of  two  synthetic  techniques;  the  Clark, 
procedure  and  the  Snyder  procedure.  The  Clark  procedure  develops  a 
unit  hydrograph  from  a  coefficient  describing  the  subbasin  time  of 
concentration  and  a  coefficient  describing  the  subbasin  natural  storage 
characteristics.  A  time-area  diagram  (sometimes  termed  a  time  delay  histo¬ 
gram  or  translation  hydrograph)  with  a  base  time  equal  to  the  time  of 
concentration  is  routed  through  a  linear  reservoir  characterized  by 
the  storage  coefficient.  The  time-area  relation  may  be  derived  from 
subbasin  physiographic  data  and  supplied  directly  or  the  program  may 
be  permitted  to  compute  and  use  an  idealized  relationship  that  consists 
of  a  simple  reverse  parabola  symmetrical  about  the  center  of  the  time 
base. 

If  a  unit  hydrograph  that  conforms  to  specified  Snyder  coefficients 
is  desired,  it  is  established  by  successive  approximations  for  the  corres¬ 
ponding  Clark  coefficients. 

Base  Flow  and  Recession.  The  program  assumes  that  the  flood  to  be 
computed  occurs  after  a  previous  flood  and  therefore  begins  at  a  flow 
on  the  recession  limb  of  the  previous  flood  hydrograph.  The  recession 
is  assumed  to  be  described  by  a.  exponential  function.  The  parameters 
required  to  describe  the  base  flow  and  recession  are  the  starting  flow, 
the  flow  at  which  recession  begins,  and  the  recession  constant.  Figure 
3  shows  the  base  flow  separation  concept  used  by  the  program. 
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Fig.  2.  General  Loss  Rate  Function  on  Snow-Free  Ground 
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Fig.  3.  Base  Flow  Separation 


Routing.  The  routing  procedures  available  In  the  program  are  the 


simplified  storage  routing  procedures  that  neglect  Inertia  effects. 

All  these  procedures  can  be  used  for  streamflow  routing  and  a  few  are 
useful  for  simple  reservoir  routing. 

The  Muskingum  procedure  determines  reach  outflow  based  on  Inflow  and 
coefficients  describing  the  reach  travel  time  and  attenuation  character¬ 
istics.  The  Tatum  method  is  a  simplification  of  the  Muskingum  procedure. 

The  modified  Puls  (sometimes  termed  storage  indication)  is  simply  a  solution 
of  the  continuity  equation  when  storage  and  outflow  are  uniquely  related 
and  is  appropriate  for  both  streamflow  and  simple  reservoirs.  The  time- 
of- storage  procedure  is  primarily  a  simple  reservoir  procedure  and  the 
working  R  and  D  is  a  combination  of  the  Muskingum  technique  and  the 
modified  Puls  procedure.  The  straddle-stagger  technique,  a  strictly 
empirical  procedure,  is  also  included. 

Automatic  Derivation  of  Parameters 
HEC-1  can  derive  loss  rate  and  unit  hydrograph  coefficients  or 
routing  coefficients  automatically.  In  order  to  optimize  coefficients, 
a  hydrograph  of  observed  runoff  must  be  supplied.  In  the  case  of  routing 
optimization,  a  pattern  hydrograph  for  intermediate  runoff  (between  given 
inflow  and  outflow)  must  also  be  supplied.  This  pattern  hydrograph  is 
automatically  multiplied  by  a  ratio  to  make  the  volume  of  intermediate 
runoff  equal  the  difference  between  inflow  and  outflow  volumes. 

Coefficients  ere  derived  by  successive  approximations  to  determine  the 
optimum  values  of  a  set  of  variables  that  will  result  in  a  minimum  value  of  an 

15 


objective  function.  The  objective  function  for  optimal  reconstitution 
is  the  weighted  root-mean-square  errors  between  computed  and  observed 
flows.  In  order  to  Improve  the  reproduction  of  peak  flows,  errors 
associated  with  high  flows  are  weighted  heavier  than  those  associated 
with  low  flows.  Each  error  square  is  multiplied  by  (Q  +  Q)/2Q  where 
Q  is  the  average  flow.  A  volume  check  is  included  in  the  hydrograph 
reconstruction  to  assure  approximate  correspondence  in  volume  between 
the  observed  and  computed  hydrographs. 

The  optimization  process  terminates  after  a  specified  number  of  steps. 

If  a  reproduction  is  not  satisfactory,  considerable  improvement  can  be  made 
in  a  second  run  by  a  routine  that  temporarily  distorts  the  observed  hydro¬ 
graph  to  force  a  better  reproduction  without  impairing  the  validity  of  the 
results.  For  example,  if  a  portion  of  a  reconstituted  hydrograph  is  too 
low,  it  can  be  fitted  better  by  increasing  a  key  flow  by  about  double  the 
discrepancy.  These  temporary  adjustments  to  the  flow  are  removed  before 
the  hydrograph  is  printed. 

Hydrograph  Balancing 

A  hydrograph  balance  routine  is  included  to  convert  any  given  hydro¬ 
graph  to  one  having  specified  volumes  within  given  durations.  Starting 
with  the  shortest  duration  specified,  the  period  of  maximum  flow  of  the 
given  hydrograph  is  determined,  and  the  sum  of  all  flows  (that  have  not 
already  been  used  in  shorter-duration  computation)  within  each  period  is 
computed.  Flows  are  adjusted  within  each  period  by  multiplying  by  the  ratio 
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required  to  obtain  the  incremental  volume  needed.  Since  the  changed 
shape  of  the  hydrograph  can  alter  the  location  of  maximum  flows,  this 
process  is  repeated  until  all  volumes  are  within  1  percent,  using 
the  derived  hydrograph  as  the  new  pattern  hydrograph  each 
iteration. 

General  Input  Structure 

HEC-1  is  designed  to  accept  input  in  card  form  that  will  describe 
the  type  of  job  to  be  performed,  its  scope  and  the  detailed  information 
required  for  processing. 

A  single  job  to  be  processed  by  HEC-1  consists  of  one  set  of  A 
through  J  cards  followed  by  repeated  sets  of  K  cards  through  Z  cards 
as  necessary.  There  are  five  types  of  jobs  that  can  be  processed  by 
HEC-1,  and  these  are  specified  by  a  key  variable. 

The  values  of  key  variables  determine  the  data  to  be  processed.  Only 
the  cards  required  for  the  desired  process  are  to  be  used,  and  they  must 
be  in  the  proper  sequence. 

Cards  K  through  Z  are  used  to  specify  the  data  relative  to  the  individual 
subarea,  routing  reach,  or  combining  operations  and  economic  analyses.  The 
order  of  computing  hydrographs  is  vital  in  certain  respects  since  the  most 
recently  computed  hydrographs  remaining  in  storage  are  the  ones  that  may  be 
combined  or  routed. 
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The  type  of  information  contained  on  the  data  cards  is  summarized  below 


Card  Code 
A 
B 
C 
D 


E,  F 
G,  H 

I 

J 

K 

L 

M-X 


Y 


Description 

Job  title 

Job  specifications 

Observed  hydrograph  to  be  reconstituted 

Routing  optimization  criteria  and  observed 
inflow  hydro graph 

Unit  graph  and  loss  rate  optimization  criteria 

Station  precipitation  data  for  all  subbasins 
of  the  watershed 

Precipitation  depth-drainage  area  data 
Multiflood,  multiplan  data 

Computation  specification  for  model  building 

Hydrograph  balancing  criteria 

Subarea  runoff  computation  data  including  pre¬ 
cipitation,  losses,  and  unit  graph  information 

Individual  reach  routing  criteria 

Economic  and  flood  frequency  data 


In  modelling  a  watershed  (nonoptimization  jobB)  the  K  cards ,  followed 
by  the  appropriate  runoff  or  routing  cards,  will  be  repeated  as  many  times 
as  necessary  in  order  to  compute  the  runoff,  routing  and  combining  of  hydro- 
graphs  in  a  logical  progression  through  the  basin. 

An  example  of  program  input  and  output  for  a  subarea  runoff  hydrograph 
computation  is  contained  on  the  following  pages. 


Id 


sub-area  Runoff  computation 
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Hardware  and  Software  Requirements 
HEC-1  has  been  developed  and  tested  primarily  on  the  UNIVAC  1108 
and  the  Control  Data  Corporation  6600  computer  systems.  It  was  then 
adapted  for  use  on  the  GE  400  series  computers.  The  program  is  written 
in  FORTRAN  IV  and  contains  about  3,000  FORTRAN  statements.  Table  1 
shows  the  hardware  requirements  and  selected  running  times  for  the  pro¬ 
gram. 


Table  1.  Hardware  Requirements  and  Running  Times 

FORTRAN  IV  Compiler 
Four  tape  and/or  disk  units* 


UNIVAC  1108 

CDC  7600 

GE  400 

Memory  (words) 

38,700 

35,400 

32,000 

Printer  (positions) 

132 

132 

120 

Compilation  (CPU  seconds) 

30 

3 

— 

Execution  of  Test  5 
(CPU  seconds) 

11 

1 

„ 

*May  only  require  two  tapes  or  disks  if  output  hydrographs  are  not  to  be 
saved  or  read  in  from  previous  jobs. 

WATER  SURFACE  PROFILES  (HEC-2) 

General  Capabilities 

HEC-2  has  been  designed  to  perform  steady  gradually  varied  flow  pro¬ 
file  computations  for  natural  channels.  The  program  accommodates  local 
obstructions  such  as  weirs,  culverts  and  bridges  so  that  a  continuous 
(uninterrupted  computation)  profile  can  be  computed  for  either  subcritical 
or  supercritical  flow.  Reference  3  contains  detailed  information  on  the 
capabilities  and  data  preparation  requirements. 
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Profile  computations  for  channel  reaches  are  based  on  solution  of  the 
one-dimensional  Bernoulli  equation  with  energy  loss  due  to  friction 
evaluated  with  the  Manning's  equation  and  other  losses  evaluated  by  appli¬ 
cation  of  the  shock-loss  equation.  The  solution  procedure  used  is 
generally  referred  to  as  the  standard-step  method. 

The  flow  field  at  a  cross  section  is  divided  into  main  channel  and 
overbank  areas  and  the  overbank  areas  are  subdivided  to  account  Tor 
nonuniform  velocity  distributions.  The  reach  distances  to  adjacent  cross- 
sections  may  be  specified  for  the  channel  and  overbank  areas. 

Energy  losses  at  obstructions  can  be  accommodated  by  three  alternative 
means.  The  loss  (rating  curve)  may  be  supplied  directly.  The  "normal" 
bridge  routine  evaluates  energy  losses  by  normal  standard-step  computations 
with  corrections  for  flow  area  and  wetted  perimeter.  The  ''special''  bridge 
routine  applies  the  principles  of  conservation  of  momentum  to  evaluate 
depths  and  thus,  indirectly,  energy  loss.  It  observes  the  hydraulic 
control  concept  of  critical  depth  in  determining  the  flow  characteristics. 
Tills  routine  also  computes  energy  losses  based  on  the  weir  equation  when 
flow  is  over  the  road  deck,  and  the  orifice  equation  when  flow  is  confined 
to  the  bridge  opening  and  under  pressure.  Bridge  characteristics  are 
described  by  the  natural  valley  cross  section,  the  top  of  roadway  and  low 
chord  profiles,  and  appropriate  discharge  and  loss  coefficients. 

The  program  has  been  carefully  designed  to  ease  data  handling  and 
manipulation  requirements.  Many  routines  and  options  are  available  for 
describing  cross  sections,  specifying  the  portion  of  the  cross  section 


Chat  Is  effective  In  passing  flow,  manipulating  the  cross  section  by 
skewing  and  raising  and  lowering  as  desired.  Utility  features  permit 
processing  many  profiles  during  a  single  run,  plotting  output,  summarizing 
computations  and  editing  data.  A  large  number  of  diagnostic  notes  are 
output  during  computations. 

The  program  has  a  few  special  routines  designed  to  assist  in 
flood  plain  management  and  flood  insurance  studies  wherein  encroach¬ 
ments  must  be  evaluated  and  floodways  must  be  designated. 

Computational  Procedures  for  Channels 

The  standard-step  method  of  computing  water  surface  elevation  at 
specific  locations  requires  a  trial  solution  of  the  Bernoulli  equation. 
Figure  4  indicates  a  schematic  of  the  computation  procedure.  The 
terms  in  the  equation  that  require  special  attention  in  natural 
channels  are  the  velocity  head  and  the  friction  loss. 

Velocity  Head.  Because  the  Bernoulli  equation  is  one-dimensional, 
the  velocity  head  term  must  be  corrected  to  account  for  the  kinetic  energy 
content  of  the  moving  total  fluid  considering  the  nonuniformity  of  the 
velocity  distribution.  The  procedure  used  is  to  subdivide  the  cross 
section  at  each  coordinate  point,  compute  the  velocity  head  for  each  sub¬ 
section  and  weight  the  incremental  velocity  head  in  accordance  with  the 
discharge  in  each  subsection.  The  velocity  for  each  subsection  is  computed 
by  the  Manning  equation  for  the  roughness  and  wetted  perimeter  of  that 
subsection. 
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Computation  Procedure 
Conditions  are  known  @  1 


•  Assume  water  surface  elevation  at@  -  WSEL  (2) 

*  ***  Compute  conveyance —  K(2)  =  IK  +  K,  +  .  . 

.  „  1.49  . _2/3  a  b 

where:  K  -  - AR 

n 

#  ***  Compute  representative  rate  of  friction  loss 

k  .  Mj i±  mi. .  s  „  2L_ 


•  Compute  WSEL(2)  „2  2 

v  -  I  z  t  V 

WSEL(2)  -  WSEL(l)  +  a  jg  +  S|v>  +  C  •  — -  a-^ 
Compare  and  repeat  until  criterion  met. 


Fig.  4.  Schematic  of  Computation  Procedure 


Friction  Loss.  The  friction  loss  between  adjacent  cross  sections  is 
computed  as  the  product  of  the  representative  rate  of  friction  loss 
(friction  slope)  and  the  weighted  reach  length.  The  weighted  reach  length 
is  determined  by  weighting  the  channel  length  for  each  overbank  area  and 
the  main  channel  by  the  discharge  flowing  in  each  of  these  major  section 
elements.  The  representative  rate  of  friction  loss  is  determined  by 
averaging  the  conveyances  of  the  adjacent  cross  section  and  then  computing 
an  equivalent  energy  grade  line  slope  from  the  Hanning  equation.  The  con¬ 
veyance  of  each  cross  section  is  the  sum  of  the  conveyances  of  each  subsection 
within  the  cross  section.  The  conveyance  of  each  subsection  is  computed 
from  the  geometric  and  hydraulic  parameters  of  the  Manning  equation  as 
indicated  in  figure  4. 

Other  Losses.  The  program  considers  other  losses  to  be  expansion  and 
contraction  losses  that  are  generated  because  of  changes  in  cross  section 
geometry.  The  loss  is  computed  as  the  product  of  the  appropriate  loss 
coefficient  and  the  absolute  value  of  the  difference  in  velocity  head. 

The  loss  is  computed  each  time  a  new  cross  section  is  encountered,  including 
those  through  bridges. 

Critical  Depth.  The  program  assumes  the  profile  being  computed  is 
either  all  subcritical  or  all  supercritical.  Either  can  be  processed  but 
must  be  done  so  separately.  Critical  depth  therefore  cannot  be  crossed 
during  a  profile  computation.  If  the  search  procedure  determines  that 
the  depth  being  computed  lies  on  the  other  side  of  critical  depth  from  the 
initial  specification,  the  program  assumes  critical  depth  at  that 
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section,  prints  a  note  and  proceeds  to  the  next  section.  Critical  depth 
is  calculated  for  every  cross  section  for  supercritical  profiles  but 
only  when  necessary  to  determine  if  crossing  is  possible  for  subcritical 
profiles.  Critical  depth  is  computed  by  a  search  procedure  that  seeks 
the  depth  that  results  in  the  minimum  value  of  specific  energy.  The 
nonuniformity  of  velocity  distribution  in  the  cross  section  is  considered 
in  th’i  computations. 

Computation  Controls.  The  solution  for  critical  depth  and  for  the 
water  surface  elevation  at  each  cross  section  requires  trial  and  error 
iterative-type  computations.  The  number  of  iterations  used  for  each 
computation  la  limited  so  that  failure  to  converge  would  not  cause  profile 
computation  to  cease.  In  those  instances  where  convergence  does  not 
occur,  messages  are  generated  to  focus  attention  on  problem  areas.  Failure 
to  converge  is  not  necessarily  caused  by  errors  since  complex  natural 
cross  sections  can  give  rise  to  Inconsistencies  in  solution  of  the  one¬ 
dimensional  equations. 

Bridge  Losses 

Energy  losses  caused  by  structures  such  as  bridges  and  culverts  are 
computed  in  two  parts.  First,  the  losses  due  to  expansion  and  contraction 
of  the  cross  section  on  the  upstream  and  downstream  sides  of  the  structure 
are  computed.  Secondly,  the  loss  through  the  structure  itself  is  computed 
by  either  the  normal  bridge  routine  or  the  special  bridge  routine. 

Normal  Bridge  Routine.  The  normal  routine  handles  the  cross  section 
at  the  bridge  just  as  it  would  any  river  cross  section  with  the  exception 
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that  the  area  of  the  bridge  below  the  water  surface  is  subtracted  from 
the  total  area  and  the  wetted  perimeter  is  increased  where  the  water 
surface  elevation  exceeds  the  low  chord.  The  bridge  deck  is  described 
by  entering  the  elevation  of  the  top  of  roadway  and  low  chord,  or  by 
specifying  a  table  of  roadway  elevation  and  station  and  corresponding  low 
chord  elevations.  When  only  top-of-roadway  and  low  chord  elevations 
are  used,  these  elevations  are  extended  horizontally  until  they  intersect  the 
ground  line.  Pier  losses  are  accounted  for  by  the  increased  wetted  peri¬ 
meter  of  the  piers.  The  normal  routine  is  particularly  applicable  for 
bridges  without  piers,  bridges  under  high  submergence,  and  for  low  flow 
through  circular  and  arch  culverts.  Whenever  flow  crosses  critical 
depth  in  a  structure,  the  special  bridge  routine  should  be  used.  The 
normal  bridge  is  automatically  used  by  the  computer,  even  though  data  was 
prepared  for  the  special  bridge  routine,  for  bridges  without  piers  and 
under  low  flow  control. 

Special  Bridge  Routine.  The  special  bridge  routine  can  be  used  for 
any  bridge,  but  should  be  used  for  trapezoidal  bridges  with  piers  where  low 
flow  controls,  for  pressure  flow  through  circular  or  arch  culverts,  and 
whenever  flow  passes  through  critical  depth  when  going  through  the  structure. 
The  special  bridge  routine  computes  losses  through  the  structure  for  low 
flow,  weir  flow  and  pressure  flow  or  for  any  confcination  of  these.  The 
type  of  flow  is  determined  by  a  series  of  comparisons  as  shown  on  figure  5 
and  as  described  below.  First,  the  energy  grade  line  elevations  are 
computed  assuming  alternately  low  flow  and  pressure  flow  control.  The 
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FIG.  5  ILLUSTRATION  OF  BRIDGE  FLOW  TYPES 


higher  energy  grade  line  elevation  determines  the  appropriate  type  of 
flow.  If  pressure  flow  appears  to  control  and  the  energy  grade  line  is 
above  the  minimum  top-of-roadway  elevation,  then  a  combination  of  pressure 
flow  and  weir  flow  exists.  If  the  energy  gradient  is  below  the  minimum 
top  of  roadway,  then  pressure  flow  alone  controls.  If  low  flow  appears  to 
control,  and  the  corresponding  energy  gradient  elevation  is  above  the  minimum 
top-of-roadway  elevation,  then  a  combination  of  low  flow  under  the  bridge 
and  weir  flow  over  the  roadway  approach  exists;  if  the  energy  elevation 
is  below  the  minimum  top-of-roadway,  then  low  flow  controls. 

Low  flow  is  further  classified  as  ClaBS  A,  B  and  C  depending  on  whether 
subcritlcal,  critical,  or  supercritical  flow  occurs  between  bridge  piers. 
Class  A  flow  (subcritlcal  flow  through  bridge)  is  solved  from  Yamell's 
energy  equation.  Class  B  and  C  flows  (critical  and  supercritical  flows, 
respectively)  are  handled  by  momentum  relations  which  equate  the  momentum 
flux  at  adjacent  locations  and  solve  for  depth.  The  flow  class  (A,  B  or  C) 
is  determined  by  comparison  of  the  momentum  flux  within  the  constriction 
based  on  upstream  and  downstream  conditions  and  the  momentum  flux  assuming 
critical  depth  in  the  constriction.  Weir  flow  is  computed  by  the  standard 
broad-crested  weir  equation.  The  approach  velocity  is  included  by  using 
the  energy  grade  line  elevation  in  lieu  of  the  upstream  water  surface 
elevation  for  computing  the  head,  H.  Where  submergence  by  tailwater  exists, 
the  discharge  coefficient  is  reduced  by  the  computer  program  according 
to  a  specified  procedure.  The  total  flow  is  computed  by  dividing  the 
weir  flow  into  subsections,  confuting  the  incremental  flow  for  each 
subsection  and  summing  all  subareas. 
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Pressure  flow  computations  use  the  orifice  flow  equation  with  the 


head  defined  as  the  difference  between  the  energy  gradient  elevation 
upstream  and  tailwater  elevation  downstream.  The  total  loss  coefficient 
K,  representing  losses  between  the  cross  sections  immediately  upstream 
and  downstream  of  the  bridge,  is  equal  to  the  sum  of  loss  coefficients 
for  intake,  intermediate  piers,  friction,  exit  and  other  minor  losses. 

Often  combinations  of  these  three  basic  types  of  flow  occur.  In 
these  cases,  a  trial  and  error  procedure  is  used  with  the  equations 
just  described  to  determine  the  amount  of  each  type  of  flow.  The  procedure 
consists  of  assuming  energy  elevations  and  computing  the  total  discharge 
until  the  computed  discharge  equals,  within  one  percent,  the  discharge 
desired. 

Program  Input  for  Profile  Computation 
A  major  portion  of  the  programming  in  HEC-2  is  devoted  to  providing 
a  large  variety  of  input  and  data  manipulation  options.  The  program 
objective  is  quite  simple — compute  water  surface  elevations  at  all  locations 
of  interest  for  given  flow  values.  The  data  needed  to  perform  those 
computations  include  flow,  starting  elevation,  cross  section  geometry  and 
roughness ,  and  reach  lengths.  The  options  available  for  providing  and 
manipulating  input  are  discussed  in  the  following  pages. 

Flow.  The  river  flow  may  be  specified  and  altered  in  several  ways. 

The  starting  flow  is  normally  specified  as  a  single  value  when  only  one 
flow  is  anticipated.  If  it  is  desired  to  use  different  flows  for  subsequent 
jobs  using  the  same  cross  sections,  a  discharge  table  may  be  used.  The  flow 
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may  also  be  changed  at  any  cross  section  such  as  at  a  confluence  with 


another  river  or  stream.  The  flow  change  at  any  cross  section  becomes 
permanent  for  subsequent  cross  sections.  The  flow  for  the  entire  profile 
may  be  Increased  or  decreased  by  a  factor  if  desired.  This  feature 
can  be  very  useful  in  sensitivity  studies. 

Starting  Elevation.  The  water  surface  elevation  for  the  beginning 
cross  section  may  be  specified  in  one  of  three  ways:  (1)  as  critical 
depth  lii  which  case  it  will  be  computed,  (2)  as  a  known  elevation,  (3)  by 
the  slope  area  method.  For  beginning  by  the  slope  area  method,  the  estimated 
slope  of  the  energy  grade  line  must  be  provided  and  the  initial  estimate  of 
the  water  surface  elevation  is  also  needed.  The  flows  computed  for  the 
fixed  slope  and  estimated  depth  are  compared  with  the  starting  flow  and 
the  initial  depth  is  adjusted  until  the  computed  flow  is  within  1  percent 
of  the  starting  flow. 

River  Flow  Geometry.  Cross  sections  are  required  at  representative 
locations  throughout  the  river  reach.  These  are  locations  where  changes 
occur  in  slope,  cross  sectional  area,  or  channel  roughness;  locations 
where  levees  begin  or  end;  and  at  bridges.  In  general,  for  rivers  of 
flat  slope  and  fairly  uniform  section  (drop  of  3  or  4  feet  per  mile) 
cross  sections  should  be  taken  at  least  every  mile.  For  steeper  slopes  and 
very  Irregular  cross  sections,  four  or  five  cross  sections  per  mile  may 
be  necessary.  Where  an  abrupt  change  occurs  in  the  cross  section,  several 
cross  sections  should  be  used  to  describe  the  change  regardless  of  the 
distance.  Every  effort  should  be  made  to  obtain  cross  sections  that 
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accurately  represent  the  river  geometry. 

Each  cross  section  in  the  reach  is  described  by  coordinates  that 
give  a  station  number  corresponding  to  the  horizontal  distance  from 
the  first  point  on  the  left  and  the  corresponding  elevation  of  the 
ground  surface.  Cross  sections  may  be  oriented  looking  either  upstream 
or  downstream  since  the  program  considers  the  left  side  to  be  the  lowest 
station  number  and  the  right  side  the  highest.  The  left  and  right 
stations  separating  the  channel  from  the  overbank  areas  must  be  specified. 
End  points  of  a  cross  section  that  are  too  low  (below  the  computed  water 
surface  elevation)  will  automatically  be  extended  vertically  by  the 
program,  if  needed,  and  a  message  giving  the  vertical  distance  extended  will 
be  printed. 

There  are  times  when  the  user  wishes  to  use  the  previous  cross 
section  as  the  current  one  (for  uniform  channels)  ,  with  or  without  a  modi¬ 
fication,  or  to  modify  the  current  cross  section  (perhaps  the  surveyed 
cross  section  is  moved  upstream  or  downstream) .  The  horizontal  dimensions 
of  the  previous  or  current  cross  section  can  be  increased  or  decreased 
by  a  factor  and  all  the  elevations  of  the  previous  or  current  cross 
sections  can  be  raised  or  lowered  by  a  constant. 

The  existing  cross  section  can  be  modified  due  to  the  excavation  of 
a  trapezoidal  channel.  The  coordinate  points  are  modified  due  to  the 
excavation,  but  no  fill  is  used.  The  bank  elevations  and  stations  are 
modified  if  the  channel  daylights  outside  the  original  bank  stations. 

If  the  alignment  of  the  excavated  channel  is  such  that  two  separate 
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channels  exist,  the  division  between  overbank  and  channel  will  be  based 
on  the  excavated  channel,  and  the  old  channel  will  be  considered  as 
overbank  (no  fill) .  It  may  be  necessary  to  change  the  reach  lengths  for 
this  case. 

Levees  require  special  consideration  in  comput’v.g  water  surface  pro¬ 
files  because  of  possible  overflow  into  areas  outside  the  main  channel. 
Informally  the  computations  are  based  on  the  assumption  that  all  area  below 
the  water  surface  elevation  is  effective  in  passing  the  discharge.  How¬ 
ever,  if  the  water  surface  elevation  is  less  than  the  top-of-levee 
elevation,  and  if  the  water  cannot  enter  the  overbanks  upstream  or  downstream 
of  that  cross  section,  then  all  flow  area  in  these  overbanks  should  not  be 
used  in  the  computations.  By  setting  a  code,  the  program  will  consider 
only  flow  confined  by  the  levees,  unless  the  water  surface  elevation  is 
above  the  top  of  one  or  both  sides  of  the  levee,  in  which  case  flow  area 
or  areas  outside  the  levee  will  be  included.  It  is  important  for  the 
user  to  study  carefully  the  flow  pattern  of  the  river  where  levees  exist 
to  determine  such  items  as  if  a  levee  were  open  at  both  ends  and  flow 
could  pass  behind  the  levee  without  overtopping  it.  Also,  assumptions 
regarding  effective  flow  areas  may  change  with  changes  in  flow  magnitude. 
Where  cross  section  elevations  outside  the  levee  are  considerably  lower 
than  the  channel  bottom,  it  may  be  necessary  to  confine  the  flow  to  the 
channel. 

Sometimes  it  is  necessary  to  insert  cross  sections  between  those 
specified  because  the  change  flow  in  characteristics  between  cross  sections 
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is  Coo  great  Co  accurately  determine  the  energy  losses.  Up  to  three 
interpolated  cross  sections  will  be  generated  between  given  cross  sections 
under  certain  conditions.  It  is  possible  to  suppress  program  generated 
interpolated  cross  sections,  and  they  should  be  omitted  when  computing 
several  profiles  on  the  same  stream  in  order  to  use  exactly  the  same 
cross  sections.  The  distance  between  cross  sections  used  in  the  computation 
can  be  specified  for  the  left  overbank,  right  overbank,  and  channel, 
respectively.  There  are  conditions  where  they  will  differ,  such  as  at 
river  bends,  or  where  the  channel  meanders  considerably  and  the  over¬ 
banks  are  straight. 

Channel  Roughness.  Since  Manning's  coefficient  of  roughness,  n, 
depends  on  such  factors  as  type  and  amount  of  vegetation,  channel 
configuration  and  stage,  several  options  are  available  to  vary  n.  The 
normal  situation  is  when  three  n  values  are  sufficient  to  describe  the 
channel  and  overbank  roughness.  Any  of  the  n  values  may  be  permanently 
changed  at  any  cross  section.  Often  three  values  are  not  enough  to 
adequately  describe  the  lateral  roughness  variation  in  the  overbanks  in 
which  case  up  to  19  n  values  and  corresponding  cross  section  stations  are 
provided.  The  n  values  thus  specified  apply  only  to  the  overbanks,  while 
the  value  used  for  the  channel  is  as  normally  entered. 

Data  Indicating  the  variation  of  Manning's  n  with  river  stage  may 
be  used  in  the  program.  This  option  applies  only  to  the  channel  area. 

It  is  possible  for  subsequent  runs  of  the  same  job  to  multiply  the 
n  values  specified  by  a  multiplier.  The  desired  multiplier  is  simply  entered 
for  each  job.  This  feature  may  be  useful  in  calibration  and  sensitivity  studies. 


Manning's  n  can  be  computed  from  known  high  water  marks  along  the 
river  reach  if  the  discharge,  relative  ratios  of  the  n  values  for  the 
channel  and  overbanks  and  the  water  surface  elevations  at  each  cross 
section  are  known.  The  "best  estimate"  of  n  for  the  first  cross  section 
must  be  provided  since  it  is  not  possible  to  compute  an  n  value  for 
this  cross  section.  The  relative  ratio  of  n  between  channel  and  over¬ 
bank  is  set  by  the  first  cross  section  and  will  be  used  for  all  subsequent 
cross  sections  unless  changed.  The  n  values  required  to  match  the  high  water 
marks  are  then  computed  provided  an  adverse  slope  in  energy  grade  line 
is  not  encountered,  in  which  case  computations  restart  using  n  values 
from  the  previous  section. 

Utility  Features  and  Output 

Several  profiles  may  be  computed  using  the  same  cross  section  data. 

A  summary  printout  will  provide  a  concise  listing  of  the  key  results  for 
all  profiles  for  each  cross  section.  The  user  may  select  a  number  of  the 
variables  to  be  retained  for  the  summary  printout. 

Multiple  Stream  Profiles.  The  water  surface  profile  computations 
may  be  extended  up  both  forks  of  a  river  or  throughout  a  whole  river  basin 
for  single  or  multiple  profiles  in  a  single  computer  run.  The  profile  is 
first  computed  for  reach  1  from  the  most  downstream  point  to  the  end  of 
one  tributary.  The  data  for  a  second  tributary  (reach  2),  whose  starting 
water  surface  elevation  was  determined  when  reach  1  was  calculated, 
follows  the  data  for  reach  1.  When  a  negative  section  number  is 
encountered,  the  program  will  search  its  memory  for  the  computed  water 
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surface  elevation  that  corresponds  to  the  negative  section  number.  It 
will  then  start  computing  the  profile  for  reach  2  with  the  previously 
determined  water  surface  elevation. 

Storage-Outflow  Data.  Punched  cards  can  be  obtained  from  HEC-2 
for  stream  routing  by  the  modified  Puls  method  using  the  program  HEC-1. 

This  option  can  be  used  only  if  multiple  profiles  are  computed  from  the 
same  cross  sectional  data.  It  may  not  be  wise  to  use  interpolated  cross 

i 

sections  since  a  different  number  of  cross  sections  might  be  interpolated  i 

between  two  given  cross  sections  for  different  magnitudes  of  discharge 

which  could  cause  inconsistencies  in  the  incremental  storage  volumes.  \ 

i 

Encroachment  Options.  Four  methods  of  specifying  encroachments  for  j 

floodway  studies  can  be  used.  Stations  and  elevations  of  the  left  and/ 

or  right  encroachment  can  be  specified  for  individual  cross  sections  as  i 

desired.  Stations  can  also  be  specified  differently  for  each  profile. 

A  fixed  topwidth  can  be  specified  which  will  be  used  for  all  cross  sections  1 

until  changed.  The  left  and  right  encroachment  stations  are  made  equi¬ 
distant  from  the  centerline  of  the  channel,  which  is  half-way  between  j 

the  left  and  right  bank  stations.  Different  topwidths  may  be  specified 

for  each  profile  of  a  multiple  profile  run.  Encroachments  can  be  specified  ! 

by  percentages  which  indicate  the  desired  proportional  reduction  in  the  natural  I 

(first  profile)  discharge  carrying  capacity  (conveyance)  of  each  cross  ! 

section.  For  example,  input  could  specify  that  for  the  second  profile  j 

I 

i 

5  percent  of  the  flow-carrying  capacity  based  on  the  first  profile, 
will  be  eliminated  on  each  side  of  the  main  channel  as  long  as 

! 
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Che  encroachments  do  not  fall  vichin  the  main  channel.  If  one  side 
cannot  carry  the  S  percent  reduction,  a  reduction  of  more  than  5  percent 
will  be  attempted  on  the  other  side.  The  first  profile  is  for  natural 
conditions;  different  sets  of  ratios  can  be  specified  for  all  subsequent 
profiles. 

Encroachments  can  be  determined  so  that  each  modified  cross  section 
will  have  the  same  discharge  carrying  capacity  (at  some  higher  elevation) 
es  the  natural  cross  section.  This  higher  elevation  is  specified  as  a 
fixed  amount  above  the  natural  (e.g.,  100-year)  profile.  The  encroach¬ 
ments  are  determined  so  that  an  equal  loss  of  conveyance  (at  the  higher 
elevation)  occurs  on  each  side  of  the  channel,  if  possible,  as  before. 

The  horizontal  distribution  of  area,  velocity  and  discharge  can  be 
requested  for  the  overbank  areas.  If  the ^number  of  subsections  carrying 
flow  in  the  overbanks  is  less  than  11,  the  distribution  using  all  sub¬ 
sections  will  be  printed.  Otherwise,  the  distribution  will  be  based 
on  subsections  that  carry  more  than  3  percent  of  the  flow. 

Cross  Section  Plot.  Plots  on  the  printer  of  any  or  all  of  the  river 
cross  sections  to  any  scale  may  be  requested.  Vertical  and  horizontal 
scales  of  the  plot  may  be  specified  constant  for  all  cross  sections  in 
the  job . 

Profile  Plot.  This  plot  includes  the  water  surface  elevation,  the 
critical  water  surface  elevation,  energy  grade  line,  channel  invert, 
left  and  right  bank  elevations,  and  the  maximum  elevation  of  the  cross 
section  for  which  hydraulic  properties  can  be  computed.  The  vertical 
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scale  of  Che  profile  may  be  determined  by  Che  user  or  by  Che  cotquiter. 
Profiles  are  plotted  automatically  for  jobs  using  more  than  five  cross 
sections,  but  may  be  suppressed  If  desired. 

Example 

The  following  pages  contain  input  for  a  sample  profile  computation 
and  normal  output,  a  summary  printout  and  a  profile  plot. 
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Hardware  and  Software  Requirements 
HEC-2  has  been  developed  and  tested  primarily  on  the  UNIVAC  1108 
and  the  Control  Data  Corporation  6600  computer  systems.  It  has  subsequently 
been  adapted  for  other  machines.  The  program  is  written  In  FORTRAN  II  and 
IV,  and  includes  about  7,000  FORTRAN  statements.  The  following  table  shows 
an  approximate  comparison  of  memory  requirements  and  execution  speeds 
for  a  number  of  computers. 


Computer 

Memory 

Requirement 

Words 

Speeds 

Compile 

Test 

CDC  6600 

63,000 

59  sec 

15  sec 

UNIVAC  1108 

52,000 

60  sec 

50  sec 

GE  437 

29,000 

1U  min 

5  min 

IBM  360/50 

70,000 

(280,000  bytes) 

14.5  min 

1.7  min 

CDC  7600 

56,000 

8  sec 

2  sec 
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